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Abstract 



Continuous time random walk models with decoupled waiting time density 
are studied. When the spatial one jump probability density belongs to the 
Levy distribution type and the total time transition is exponential a generalized 
superdiffusive regime is established. This is verified by showing that the square 
width of the probability distribution (appropriately defined)grows as t 2 ^ with 
< 7 < 2 when t — ► oo. An important connection of our results and those of 
Tsallis' nonextensive statistics is shown. The normalized q-expectation value of 
x 2 calculated with the corresponding probability distribution behaves exactly 
as t 2 / 7 in the asymptotic limit. 
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Recently there has been an increasing interest in dynamical processes that display anoma- 
lous diffusiorJlHl These phenomena have been characterized by a non linear time dependence 
of the mean square displacement of the walker 

(x 2 (t)} ~ e (i) 

with a ^ 1 (a — 1 gives normal diffusion) since (x 2 ) is an usual estimator of the square 
width of the probability distribution at time t. In this way, for anomalous diffusion we have 
that the probability distribution width grows faster (slower) for a > 1 (a < 1) than it does 
for normal diffusion. Examples of Eq.(l) with a < 1 (dispersive diffusion or subdiffusion) 
are found in disordered media like glasses and fractals structures3~0. Chaotic dynamics 
and turbulence, on the other hand, give rise to enhanced diffusion (superdiffusion), i.e., 
ME. Diffusion like behaviour can be modeled in the framework of the continuous 
time random walk (CTRW), where the central magnitude is the waiting time density (WTD), 
9{z,t), as was first proposed by Montroll and WeissR In this way dxdt represents 

the probability that the walker makes a jump of length x after a time t since the last jump. 

When it exists the walker mean square displacement from the origin, at time t, can be 
calculated from the Fourier transform of the spatial probability density P (x, t) according to 

<x 2 (t)>= J x 2 P(x,t) dx = - 9 g^'\ =o (2) 

where P(x, t) dx is the conditional probability of finding the walker between x and x + dx 
at time t given that it started at x = at time t = 0, and P(k,t) represents its Fourier 
transform in the space variable. As was shown in previous papers, it is imperative to use a 
coupled WTD in order to have a finite step mean square displacement in the super diffusive 
regime i. e. when the distribution width grows faster than in normal diffusion. 
This requirement can be understood since for a decoupled WTD 

V(x,t) = p(x) ij)(t) (3) 

a divergent mean square displacement for the one step distribution will make infinite (x 2 (t)) 
also. This would suggest a formally infinite diffusion coefficient as pointed out by Ott et a/.0, 



while the experimental results for elongated micelles analyzed in their paper let them define 
"a (distance-dependent) diffusion constant. Nevertheless the coupled WTD requirement is 
not necessary if we consider alternative estimators of the width fo the probability distribution 
when (x 2 (t)) diverges, as we will do in this letter. 

In particular it is well known that the Levy family of probability densities^, has divergent 
second moment. In spite of this divergence of their second moment a finite width can still 
be defined for these distributions by considering alternative estimators, such as the inverse of 
P (x = 0, t) or the width at half of the maximun value. Notice that with any of these criteria 
we get a linear law for the time dependence of the square width in normal diffusion as we 
do with (x 2 (t)). One of the most familiar examples of Levy distributions is the Lorentz 
or Cauchy distribution. This distribution arises naturally in the study of the line shape of 
fluctuation spectra of stochastic Markov processes!!. It is known that the Fourier transform 
of Levy distributions are, by definition, given by0 

p(k) = exp(-c|fc| 7 ) , < 7 < 2 (4) 

In the framework of a non-extensive thermostatistics, Tsallis and collaborators^! have 
obtained, from a variational principle, probability densities that share with the Levy family 
of density functions the asymptotic behavior for large steps. The densities obtained in these 
papers are (with appropriate simple constraints) the extremal of a functional entropyi! S q 
where the parameter q < 3 identifies a particular density. Tsallis' formulation generalizes 
the classical Boltzmann-Gibbs (obtained for the particular value q = 1) variational formu- 
lation. The relationship between parameter q and the exponent 7 are related by 7 = |5f 
in the range | < q < 3 , where the second moment diverges, while 7 = 2 for q < | . 
They also showBi that the N-jump process constructed from these densities scales as N 2 ^' 
for N — > 00. 

Our purpose here is to extend the Tsallis' results to a CTRW with a decoupled WTD. 
One important result of this extension is that in case that the one step probability p(x) (|3|) 
is of the Levy type (f|), the width of the conditional probability density P(x, t) grows as t 1 ^ 1 



for t — > oo, so that for 7 < 2 we obtain a superdiffusive regime. 

Let us start by re-obtaining Tsallis' results for the N-jump distribution for a Levy flight. 
The probability density for the walker position after iV steps is given by the N-fold convo- 
lution of the one-jump distribution 

P(x, N) = p(x) * p(x) * * p(x) (5) 

whose Fourier transform is 

P(k,N) = [p(k)] N = exp(-ciV|A;| 7 ) (6) 

by making use of Eq. Anti-transforming the last equation we finally obtain 

1 r°° 

P(x,N) = — dkexp(—cN\k\' y ) exp(ikx) 

Z7T J-00 

(7) 

= N~ lh p(N~ lh x) 
where a change of variable has been made. 

Eq. (0) is valid for all N and shows that the width of the N-jump distribution, when 
appropriately defined, is that of the one-jump distribution times the factor N 1 ^. The 
prefactor ensures the normalization condition of the probability density. As it is readily 
verified, estimators such as the inverse of P (x = 0, t) or the width at half of the maximun 
value, which are finite for these distributions, exhibit the N 1 ^ growth. In the CTRW scheme 
a similar result may be obtained when the transition time probability density is 

ip(t) = Aexp(-At) (8) 

In this case 

1 r°° 

P(x,t) = — dk exp[-At (1 -p{k))] exp(ikx) (9) 
2% J-00 

satisfies the Master equation: 
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dP{x,t) 



/oo 
dx' p(x - x') P(x', t) -XP(x,t) (10) 
-oo 



dt 

and may be expressed as a power series in time as follows: 

oo /\ j.\n 

P(x,t) = exp(-Xt)5(x) + exp(-At) ]T [ — r~~ p(n~ 1/7 x)} (11) 



n=l W! 



where 5(x) is Dirac's delta. 

In the limit A t ^> 1 only small values of k contribute to the integral in Eq. (^)0, giving 
rise to 

1 r°° 

P(x,t) ~ — / dk exp(-A t c|£;| 7 ) exp(ikx) for At> 1 (12) 

By making a change of variables we get 

P{x,t) ~ {\t)- lh p[{\t)- lh x] (13) 

then the width of the probability distribution at time t grows as (A t) 1 ^ 1 in a similar way as 
was obtained for the discrete time case. As we have mentioned before in normal diffusion 
< x 2 (t) > ~ t, i.e., the square width grows as t, while, in the N-jump Levy process, as 
we have shown, the square width behaves as N 2 ^ 1 , and as t 2 ^ in the continuous time case. 
Since < 7 < 2 we conclude that a superdiffusive regime is established even in the separable 
CTRW case when the one-jump space distribution belongs to the Levy family. 

In Figure 1 we show a plot of the width of the probability distribution versus t (in 
logarithmic scale) for various values of the parameter 7. Here two criteria have been used 
to define the width of the distributions; one of them just consists in taking the inverse of 
P(x = 0,t), and the other one is the width taken at half of the maximun value of the 
distribution. We observe that, for \t > 10, the width grows linearly with t 1 / 7 , according to 
our prediction. 

It is worthwhile to notice that our results may be connected with those of Tsallis and 
collaborators obtained in the framework of nonextensive statisticsB'i. Followig these authors, 
we calculate the normalized q-expectation value of x 2 for a N-jump process 



< x 2 > q (iV) 



JZodxx 2 [P(x,N)]« 
JZodx [P{x,N)]< 



(14) 



with the following constraint imposed on the single step distribution 



a 



2 



f^o dx x 2 p(x) q 
IZodxp(x)i 



(15) 



Using Eq. (0) and performing a change of variable we obtain 



< x 2 > q (iV) = N 2h a 2 



(16) 



This result may be extended to the continuous time case for A t ^> 1 



< x 2 >q (t) = (A t) 2h a 2 



(17) 



The last two results suggest that the normalized q-expectation values of x 2 can be used as 
an estimator of the square width of the Levy type distributions obtained in!, whose normal 
(q = 1) expectation value of x 2 (second moment) diverges. 

To conclude, let us emphasize the main result of this letter: a CTRW, with a separable 
waiting time density *f?(x,t) = p(x) ip(t), whose one-jump space distribution p(x) belongs 
to a Levy family, gives rise to a generalized superdiffusive process in the case of an expo- 
nential form, Eq. @, for the total transition probability density if)(t). If, however, ip(t) has 
a long time tail the above conclusion may not be true; work along this line is in progress. 
Finally, the non-extensive thermo statistics seems to be the appropriate framework to tackle 
anomalous diffusion. 
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Figure caption 

Figure 1: Width of the probability distribution P(x,t) defined in two ways: Wi = width 
at the half height of P(x,t) and Wi —l/[P(x = 0,t] for various one-step Levy distributions 
determined by the values of 7 (see text). Notice the linear behavior of the width, with both 
criteria, in the asymptotic limit \t — > 00. The slopes of the curves are consistent with the 
value I/7 predicted by the theory. 
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